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We present a phase-plane analysis of cosmologies containing a barotropic fluid with equation of 
state P7 = (7 — l)P7i plus a scalar field </> with an exponential potential V oc exp(— A«;<^) where 
= SttG. In addition to the well-known inflationary solutions for < 2, there exist scaling 
solutions when > 87 in which the scalar field energy density tracks that of the barotropic fluid 
(which for example might be radiation or dust). We show that the scaling solutions are the unique 
late-time attractors whenever they exist. The fiuid-dominated solutions, where V{(j>)/ p-y — > at late 
times, are always unstable (except for the cosmological constant case 7 = 0). The relative energy 
density of the fluid and scalar fleld depends on the steepness of the exponential potential, which is 
constrained by nucleosynthesis to A^ > 20. We show that standard inflation models are unable to 
solve this 'relic density' problem. 
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I. INTRODUCTION 

Scalar fields have come to play a central role in cur- 
rent models of the early universe. The self-interaction 
potential energy density of such a field is undiluted by 
the expansion of the universe and hence can act like an 
effective cosmological constant driving a period of infla- 
tion. The detailed evolution is dependent upon the spe- 
cific form of the potential I^ as a function of the scalar 
field's expectation value (j). 

A common functional form for the self-interaction po- 
tential is an exponential dependence upon the scalar field. 
It is to be found in higher-order |l[ or higher-dimensional 
gravity theories Q . In string or Kaluza-Klein type mod- 
els the moduli fields associated with the geometry of 
the extra dimensions may have effective exponential po- 
tentials due to curvature of the internal spaces, or the 
interaction of moduli with form fields on the internal 
spaces. Exponential potentials can also arise due to non- 
perturbative effects such as gaugino condensation Q . 

The possible cosmological roles of exponential poten- 
tials have been investigated before, but almost always as 
a means of driving a period of cosmological inflation [^,^ . 
This requires potentials that are much flatter than those 
usually found in particle physics models. The purpose of 
this paper is to emphasize that scalar fields with expo- 
nential potentials may still have important cosmological 
consequences even if they are too steep to drive a period 
of inflation . We will present a phase-plane analy- 
sis to show that scalar fields with exponential potentials 
contribute a non-negligible energy density at nucleosyn- 
thesis unless they are unusually steep. This 'relic density' 



problem is not alleviated by standard models of inflation. 



II. AUTONOMOUS PHASE-PLANE 

We will consider a scalar field with an exponential po- 
tential energy density V — Vq exp(— Ak(/)) evolving in a 
spatially-flat Friedmann-Robertson-Walker (FRW) uni- 
verse containing a fluid with barotropic equation of state 
Pj = i"f — 1)P7, where 7 is a constant, < 7 < 2, such 
as radiation (7 — 4/3) or dust (7 = 1). The evolution 
equations for a spatially-flat FRW model with Hubble 
parameter H are 
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subject to the Friedmann constraint 
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where k? = SttG. The total energy density of a homoge- 
neous scalar field is — ^^/2 + y(0). 
We define 
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The evolution equations can then be written as a plane- 
autonomous system: 
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TABLE I. The properties of the critical points. 



a,' = _3a, + A,/^y2^^x[2x2+7(l-a:'-y')] , (6) 
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where a prime denotes a derivative with respect to the 
logarithm of the scale factor, N = ln(a), and the con- 
straint equation becomes 
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Note that from the constraint equation we have 
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This is bounded, < + < 1, tor a non- negative 
fluid density, p-y > 0, and so the evolution of this sys- 
tem is completely described by trajectories within the 
unit disc. The lower half-disc, y < 0, corresponds to con- 
tracting universes. As the system is symmetric under the 
reflection {x,y) — > (x,—y) and time reversal t — > —t, we 
only consider the upper half-disc, j/ > in the following 
discussion. 

The effective equation of state for the scalar field at 
any point is given by 
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Fixed points at finite values of x and y in the phase- 
plane correspond to solutions where the scalar field has 
a barotropic equation of state and the scale factor of the 
universe evolves as a oc where p = 2/370. 

Depending on the values of 7 and A, we have up to 
five fixed points (critical points) where x' = and y' = 
which are listed in Table I. A full analysis of the stability 
is given in the Appendix. 

Two of the fixed points {x = ±1, y — 0) correspond to 
solutions where the constraint Eq. (^ is dominated by 
the kinetic energy of the scalar field with a stiff equation 



of state, 70 — 2. As expected these solutions are unstable 
and are only expected to be relevant at early times. 

More surprisingly, however, we find that the barotropic 
fluid dominated solution {x — 0, y = 0) where fi^ = 
is unstable for all values of 7 > 0. We will discuss the 
critical case where 7 = later. But for any 7 > 0, and 
however steep the potential (i.e. whatever the value of 
A), the energy density of the scalar field never vanishes 
with respect to the other matter in the universe. 

We are left with only two possible late-time attractor 
solutions. One of these is the well-known scalar field 
dominated solution (fi^ = 1) which exists for sufficiently 
flat potentials, A^ < 6. The scalar field has an effective 
barotropic index 7^ = /3 giving rise to a power-law 
inflationary expansion (a > 0) for A^ < 2. Previous 
phase-plane analyses 1^ have shown that a wide class of 
homogeneous vacuum models approach the spatially-flat 
FRW model for A^ < 2. We have shown that this scalar 
field dominated solution is a late-time attractor in the 
presence of a barotropic fluid when we have A^ < 37. 

However for A^ > 37 we find a different late-time at- 
tractor where neither the scalar-field nor the barotropic 
fluid entirely dominates the evolution. Instead we have 
a scaling solution where the energy density of the scalar 
fleld remains proportional to that of the barotropic fluid 
with ^6 — 37/A2. This solution was first found by Wet- 
terich jg] and shown to be the global attractor solution 
for A^ > 37 in Ref. §. 

The regions of (7, A) parameter space leading to differ- 
ent qualitative evolution are indicated in Fig. 

I. A2 < 37. See Fig. | 

Both kinetic-dominated solutions are unstable 
nodes. The fluid-dominated solution is a saddle 
point. The scalar field dominated solution is the 
late-time attractor, and is inflationary in parame- 
ter region la and non-inflationary in region lb. 

II. 37 < A2 < 6. See Fig. | 

Both kinetic-dominated solutions are unstable 
nodes. The fluid-dominated solution is a sad- 
dle point. The scalar field dominated solution is 
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FIG. 1. Regions of (7, A) parameter space, as identified in 
the text. Solutions to the left of the dotted line are inflation- 
ary. 

a saddle point. The scaling solution is a stable 
node/spiral. 

III. 6 < A2. See Fig. |. 

The kinetic-dominated solution with Aa; < is 
an unstable node. The kinetic-dominated solution 
with Aa; > is a saddle point. The fluid-dominated 
solution is a saddle point. The scaling solution is a 
stable spiral. 



III. COSMOLOGICAL CONSEQUENCES 

The possible role of scalar fields with exponential po- 
tentials with A^ < 2 in driving an inflationary expansion 
of the early universe has already received considerable 
attention and so we will not devote much time to dis- 
cussing this here. However we note that because the 
scalar field dominated solution is the late-time attractor 
for A^ < 87, the existence of such a scalar field today is 
ruled out unless its energy density has been greatly sup- 
pressed relative to the attractor value (Sl^ — 1) for most 
of the "dust-dominated" era where 7 = 1. Nonetheless 
such models have been considered as possible 'decaying 
cosmological constant' models P,p|,p^. 

The peculiar properties of the scaling solution, which 
we see is the late-time attractor for exponential poten- 
tials with A^ > 87, are sufficie ntly novel to merit greater 
investigation. (See also Refs. [[ll|J^.) 

The most striking possibility is that a scalar field 
with an exponential potential could comprise a signifi- 
cant fraction of the energy density of our universe today. 
Because the effective barotropic index of the homoge- 
neous scalar field would mimic pressureless dust, its dy- 
namical effect would be exactly like cold dark matter. 
For instance, if A = 8 then we expect fi^ = 1/8 today. 
However the inhomogeneous field can have a different 
equation of state modifying the evolution of large-scale 
structure in the universe, as has recently been investi- 
gated elsewhere |^,|lo|. 
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FIG. 2. The phase plane for 7 = 1, A = 1. The 
late-time attractor is the scalar field dominated solution with 
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FIG. 3. The phase plane for 7 = 1, A = 2. The s calar 
field dominated solution is a saddle point at a: = y^2/3, 
y = -^1/3, and the late-time attractor is the scaling solu- 
tion with X — y — -^3/8. 
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FIG. 4. The phase plane for 7 = 1, A = 3. The late-time 
attractor is the scaling solution with x — y = -y/l/G. 



The main problem with this scenario is if the scalar 
field has a significant contribution to the energy density, 
throughout the present dust-dominated era, then, unlike 
conventional cold dark matter, it should also have had 
a significant effect during the radiation-dominated era. 
For A = 8 we expect fi^ = 4/9 when = p^/S. 

The tightest constraint on the total energy density 
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of the universe comes from models of nucleosynthesis 
1^,0. The primordial abundances of the light elements 
place tight constraints on the expansion rate, and hence 
the energy density, at the time of nucleosynthesis, when 
T - 1 MeV. If we require < f^^'^'^ at the time of 
nucleosynthesis, then this implies 



(11) 



The current upper bound on 17^ at nucleosynthesis is es- 
timated to be in the range 0.13 to 0.2 [^; we'll adopt 
the higher value to be conservative. Satisfying the nucle- 
osynthesis bound requires > 20. 

Thus there is a relic abundance problem for any par- 
ticle physics theories that predict the existence of scalar 
fields with exponential potentials with A^ < 20 at low 
energies (T ~ 1 MeV) . Scalar fields with exponential po- 
tentials completely dominate the energy density of the 
universe at nucleosynthesis for A^ < 4, and still have an 
unacceptably high energy density at nucleosynthesis for 
4 < A^ < 4/^7™'''^ unless the initial energy density in the 
field is extraordinarily low. 

To quantify exactly how small the initial energy den- 
sity must be to evade the bound Eq. (pT|), we expand to 
first-order about the fluid-dominated solution to find the 
rate at which fi^ grows away from zero — see Appendix. 
We find flij, oc a^'* , which implies that the scalar field 
energy density remains essentially constant due to the 
large friction term in the evolution equation Eq. ^ as 
the barotropic fluid density redshifts as oc a~^'^ . Thus 
the scalar field acts like a cosmological constant until fi^ 
approaches its attractor value. To have not reached the 
attractor by some time tf requires an initial value at U 
satisfying 



< 



(12) 



IV. THE ROLE OF INFLATION 

The usual cosmological solution to relic abundance 
problems is a period of inflation, during which the un- 
wanted relics have their energy density redshifted to a 
negligible value relative to the potential energy U of the 
inflaton field d.^ If inflation ends at some energy density 
Piih) ~ M'^, then from Eq. ( [T^ ) we require 



^ , , . ,'1 MeV\' 
^ ( ^ 



(13) 



*Note that we are assuming that the field a is unrelated to 
the field which we have been discussing up until now. 



for the scalar fleld not to have reached the scaling at- 
tractor solution by the time of nucleosynthesis. However 
we will now show that the expected density of the scalar 
field at nucleosynthesis is not significantly affected by 
standard models of inflation. 

During conventional slow-roll inflation the infiaton 
field (T has an effective equation of state 7 « 2e/3, where 
e is the (non-negative) slow-roll parameter controlling the 
slope of the potential jl^] . We can treat 7 as effectively 
constant provided the approach to the scaling attractor 
(determined by the largest eigenvalue for linear pertur- 
bations) is faster than the movement of that attractor. 
During slow-roll inflation this requires 



< 2e 



(14) 



which is valid for most models, including chaotic inflation 
with U{a) oc ct" for n>2. 

The infiaton-dominated solution where U,^ = is not 
an attractor solution unless 7 = and so for A^ > 2e the 
attractor solution is the scaling solution with 



(15) 



Unlike ordinary matter the energy density of a scalar 
field with an exponential potential does not vanish rela- 
tive to the inflaton energy density during inflation, even 
though for A^ > 2 the exponential potential would not 
have an inflationary equation of state in the absence of 
the infiaton. 

To evade the nucleosynthesis bound on A and satisfy 
Eq. ( p^ ) requires an extremely small value of e. In slow- 
roll infiation e must be smaller than unity, but is not 
usually very small. For instance in chaotic inflation with 
a potential U{a) = r7i^(T^/2, e ~ 0.01 when perturba- 
tions on the current horizon scale were generated, and 
increases to unity at the end of inflation before the infla- 
ton starts oscillating about the minimum of its potential. 
Thus the energy density of a scalar field with an exponen- 
tial potential is not significantly reduced during chaotic 
inflation, and can rapidly attain its attractor solution 
when the radiation-dominated era begins. 

Models such as hybrid inflation [|l^ or thermal infla- 
tion are distinctive in that e can remain very small 
during the flnal stages of inflation and thus fl^f, may be 
very small at the attractor scaling solution. However 
in this case the condition Eq. (|l4| ) for quasi-constant 
7 is violated and the attractor value for fl^ approaches 
zero faster than the actual solution. We effectively have 
7^0, which corresponds to the critical case where 
the scaling solution tends to the fluid-dominated solu- 
tion with X = y = 0. 

For 7 = the largest eigenvalue for linear pertur- 
bations vanishes (see Appendix) and we must consider 
higher-order perturbations about the critical point to de- 
termine its stability. We find that a; = y = is a stable 
attractor, but that trajectories only approach this as the 
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logarithm of the scale factor, N. The late-time evolution 
is given by 



2 V6 

y X 
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Thus even the extreme case of a cosmological constant (or 
constant false-vacuum energy density) only dilutes the 
energy density of the scalar field as the logarithm of the 
scale factor, fl^ cx Thus a model such as thermal 

inflation, which is so effective at diluting the abundance 
of relic moduli particles has a negligible effect on 
the relic density of the scalar field (j) as it only lasts for a 
small number of e-foldings (typically about 15) and the 
scalar field can rapidly return to its scaling solution after 
inflation. Even in the case of hybrid inflation the relic 
density after inflation may be significant unless we have 
an extremely large number of e-foldings during inflation. 
To satisfy Eq. (O) requires 
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A2 V 1 MeV . 
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which could be of order 10^" e-foldings for a typical hy- 
brid inflation model. Requiring so much expansion would 
be a significant constraint on the model. 



V. CONCLUSIONS 



APPENDIX: STABILITY OF THE CRITICAL 
POINTS 

In order to study the stability of the critical points 
(xcUc) we expand about these points 



X ^ Xc + u , y = Uc + v , 



(Al) 



which when substituted into Eqs. (||) and (^ yields, to 
first-order in the perturbations, equations of motion 



M 



(A2) 



The general solution for the evolution of linear perturba- 
tions can be written as 

u = u+ exp(m+iV) + u_ exp(m_ A^) , (A-3) 
V = exp(m+A^) + exp(m_Af) , (A4) 

where m± are the eigenvalues of the matrix Ai . Thus for 
stability we require the real part of both eigenvalues to 
be negative. 

For the critical points listed in Table 1 we find: 

Fluid- dominated solution: 



We have presented a phase-plane analysis of the evo- 
lution of a spatially-fiat FRW universe containing a 
barotropic fluid plus a scalar field with an exponential 
potential V{(j)) = Vb exp(— Ak(/>). We have shown that the 
energy density of the scalar field dominates at late times 
for A^ < 87. For A^ > 87 we find that the barotropic 
fiuid does not completely dominate and the energy den- 
sity of the scalar field remains a fixed fraction of the total 
density at late times. 

This leads to a relic density problem at nucleosynthesis 
in such models if A^ ^ 20 . Standard models of inflation 
do not significantly dilute the initial density of the expo- 
nential potential and do not alleviate this bound. Only 
inflation models with effectively constant energy density 
and an exponentially large number of e-foldings (such as 
some models of hybrid inflation) would be able to weaken 
this bound, unless the radiation-dominated era only be- 
gins shortly before nucleosynthesis, e.g. Ref. |l^. 

We emphasize that we have assumed that there is no 
direct coupling between the exponential potential and 
other matter. The only interaction is gravitational. 



Kinetic- dominated solutions^ (x^ — ±1, yc — Q)'- 

m- = y|(V6TA) , m+ = 3(2-7) . (A6) 



Scalar field dominated solution: 
A2-6 



m,+ = A — 87 



(A7) 



Scaling solution: 

8(2-7) 



m± = 



i± Wi- 



87(A^ - 37) 
A2(2-7) 



(A8) 
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